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The experimental realization of various spin ladder systems has prompted their detailed 
theoretical investigations. Here we study the evolution of ground-state magnetization 
with an external magnetic field for two different antiferromagnetic systems: a three-legged 
spin- 1/2 ladder, and a two-legged spin- 1/2 ladder with an additional diagonal interaction. 
The finite system density-matrix renormalization-group method is employed for numerical 
studies of the three-chain system, and an effective low-energy Hamiltonian is used in the 
limit of strong interchain coupling to study the two- and three-chain systems. The three- 
chain system has a magnetization plateau at one-third of the saturation magnetization. 
The two-chain system has a plateau at zero magnetization due to a gap above the singlet 
ground state. It also has a plateau at half of the saturation magnetization for a certain 
range of values of the couplings. We study the regions of transitions between plateaus 
numerically and analytically, and find that they are described, at first order in a strong- 
coupling expansion, by an XX Z spin-1/2 chain in a magnetic field; the second-order terms 
give corrections to the XXZ model. We also study numerically some low-temperature 
properties of the three-chain system, such as the magnetization, magnetic susceptibility 
and specific heat. 



PACS number: 75.10.Jm, 75.60.Ej 



I. INTRODUCTION 



One-dimensional and quasi-one-dimensional quantum spin systems have been studied extensively in 
recent years for several reasons. Many such systems have been realized experimentally, and a variety of 
theoretical techniques, both analytical and numerical, are available to study the relevant models. Due to 
large quantum fluctuations in low dimensions, such systems often have unusual properties such as a gap 
between a singlet ground state and the excited nonsinglet states; this leads to a magnetic susceptibility 
which vanishes exponentially at low temperatures. Perhaps the most famous example of this is the 
Haldane gap which was predicted theoretically in integer spin Heisenberg antiferromagnetic chainsu, and 
then observed experimentally in a spin-1 system Ni{C2HsN2)2N02{Cl04)n. Other examples include the 
spin Ladder systems in which a small number of one-dimensional spin- 1/2 chains interact amongst each 
otheiD. It has been observed that if the number of chains is even, i.e., if each rung of the ladder (which 
is the unit cell for the system) contains an even number of spin-1/2 sites, then the system effectively 
behaves like an integer spin chain with a £ap in the low-energy spectrum. Some two-chain ladders which 
show a gap are (yO)2P20Ya, SrCu203U and Cu2{C5Hi2N2)2ClM- Conversely, a three-chain ladder 
which effectively behaves like a half-odd-integer spin chain and does not exhibit a gap is Sr2Cu305U. 
A related observation is that the quagi-one-dimensional system CuGeOs spontaneously dimerizes below 
a spin-Peierls transition temperaturdZI; then the unit cell contains two spin-1/2 sites and the system is 
gapped. 

The results for gaps quoted above are all in the absencCj-pf an external magnetic field. The situation 
becomes more interesting in the presence of a magnetic fieldu. Then it is possible for an integer spin chain 
to be gapless ancLci-half-odd-integer spin chain to show a gap above the ground state for appropriate 
values of the fieldETO. This has been demonstrated in several models using a vpJety of methods such 
as exact diagonalization of small systems, bosonization and conformal field theoryEjO, and perturbation 
theoryE3. In particular, it has been shown that the magnetization of the system can exhibit plateaus at 
certain nonzero values for some finite ranges of the magnetic field. Further, for a Hamiltonian which is 
invariant under translation by one unit cell, the value of the magnetization per unit cell is quantized to 
be a rational number at each plateau. 

The necessary (but not sufficient) condition for the magnetization quantization is given as follows§l. 
Let us assume that the magnetic field points along the z axis, the total Hamiltonian H is invariant under 
spin rotations about that axis, and the maximum possible spin in each unit cell of the Hamiltonian is 
given by S. Consider a state ip such that the expectation value of Sz per unit cell is equal to uis in that 
state, and ^l) has a period n, i.e., it is invariant only under translation by a number of unit cells equal to 
n or a multiple of n. (It is clear that if n > 2, then there must be n such states with the same energy, 
since H is invariant under a translation by one unit cell). Then the quantization condition says that a 
magnetic plateau is possible at the state ip, i.e., there is a range of values of the external field for which 
ip is the ground state and is separated by a finite gap from states with slightly higher or lower values of 
total Sz, only if 

n { S — iris ) — an integer. (1) 



This condition is very useful because it enables us to restrict our attention to some particular values of 
TTis and n when searching for possible plateaus in a given model. Note that the saturated state in which 
all spins point along the magnetic field trivially satisfies (|l|) since it has — S (or —S) and n = 1. 

In this paper, we will study the magnetization as a function of the applied field for a two- and three- 
chain laclder.. We will do so both numerically, using the density-matrix renormafLizjation group method 
(DMRG)c§E3, and perturbatively, using a low-energy effective Hamiltonian (LEH)EjS. Our analysis will 
extend the currently known results in many ways. We have used DMRG to study two-spin correlation 
functions in the ground state, and some finite-temperature thermodynamic properties such as magnetic 
susceptibility and specific heat. Further, our LEH goes up to the second order in a strong-coupling ex- 
pansion. Whenever possible, we will use the analytical results from the LEH to understand the numerical 
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results. The fkS|Uprder LEH will turn out to be the well-studied XX Z spin-1/2 chain in a longitudinal 
magnetic fiplHPtP^ and it will usually prove to be sufhcient for a qualitative understanding of the results. 
However, we will find it necessary to invoke the second-order results (which give corrections to the XXZ 
model) for a more accurate comparison with the numerics. 

The paper is organized as follows. In Sec. II, we will present all the numerical results we have obtained 
for the three-chain ladder using DMRG. We will see that there is a finite energy gap and exponentially 
decaying spin correlations at each plateau, while there is no gap and the two-spin correlations decay as 
powers in between two plateaus. We will also study how the plateaus gradually disappear and how the 
susceptibility and specific heat evolve as we increase the temperature. In Sec. Ill, we will derive the LEH 
for the same model and show how it can be used to understand some of the numerical results in Sec. II. 
We will also deiive the LEH for a two-chain ladder which can be thought of as a dimerized and frustrated 
spin-1/2 chainEH, and we will use it to understand magnetization plateaus in the ground state. We will 
see that for certain values of the dimerization and frustration, the ground state can spontaneously break 
translation invariance leading to an additional plateau at an intermediate value of the magnetization. In 
Sec. IV, we will summarize our results and point out some directions for future studies. 



II. DENSITY-MATRIX RENORMALIZATION-GROUP STUDY OF THE THREE-CHAIN 

LADDER 

We have numerically studied a three-chain spin-1/2 ladder governed by the Hamiltonian 



where a denotes the chain index, n denotes the rung index, h denotes the magnetic field (we have 
absorbed the gyromagnetic ratio g and the Bohr magneton /is in the definition of h), and J, J' > 0. 
For convenience, we choose h > since the region /i < can be deduced from it by reflection about 
the zero field. It is convenient to scale out the parameter J, and quote all results in terms of the two 
dimensionless quantities J' / J and h/J. If the length of each chain is L, the total number of sites is 
TV = 3L. Since the total is a good quantum number, it is more convenient to do the numerical 
computations without including the magnetic-field term in (^, and then to add the effect of the field at 
the end of the computation. The label n-\-\ (or a-\-\) is appropriately interpreted for periodic boundary 
conditions along the chain (or rung). 

For the ground state properties, we have only considered an open boundary condition (OBC) in the 
rung direction, namely, the summation over a in the first term of (H) runs over 1,2. However, for low- 
temperature properties, we have studied both OBC, as well as a periodic boundary condition (PBC) in 
the rung direction in which we sum over a = 1, 2, 3 in the first term. (Only the OBC is realized along 
the rungs in the experimental systems studied so far. However, PBC along the rungs is interesting for 
theoretical reasons as we will see below). 

For small systems, we have performed exact diagonalization with periodic boundary conditions in 
the chaiUj-direction. For larger systems, we have done DMRG calculations (using the finite system 
algorithmtj) with open boundary conditions in the chain direction. For exact diagonalization, we have 
gone up to 24 sites, i.e., a chain length of 8. With DMRG, we have gone up to 120 sites (chain length 
of 40) after checking that the DMRG and exact results match for 24 sites. The number of dominant 
density matrix eigenstates, corresponding to the m largest eigenvalues of the density matrix, that we 
retained at each DMRG iteration was m = 80. In fact, we varied the value of m from 60 to 100, and 
found that m = 80 gives satisfactory results in terms of agreement with exact diagonalization for small 
systems and good numerical convergence for large systems. For inputting the values of the couplings into 



3 



3 





3 



the numerical programmes, it is more convenient to think of the system as a single chain (rather than as 
three chains) with the Hamiltonian 

H ^ ^J' [1 ~ cos {^)]S,- S,+i + J ^ S, . S,+3 . (3) 

i i 

The system is grown by adding two new sites at each iteration. Note that our method of construction 
ensures that we obtain the three-chain ladder structure after every third iteration when the total number 
of sites becomes a multiple of 6. At various system sizes, starting from 48 sites and going up to 120 sites 
in multiples of 6 sites, we computed the energies after doing three finite system iterations; we found that 
the energy converges very well after three iterations. The energy data is used in Figs. ^ and |^ below. 
After reaching 120 sites, we computed the spin densities and correlations after doing three finite system 
iterations. This data is used in Figs. § - |[ 

All our numerical results quoted below are for J/ J' — 1/3. We chose this particular value of the 
ratio for two reasons; there is a particularly broad magnetization plateau at rris = 1/2 which can be 
easily found numerically, and that value of the ratio is sufficiently deep inside the strong-coupling regime 
that the second-order perturbation expansion of Sec. II gives results which compare very well with the 
numerics. 

We now describe the various ground-state properties we have found with OBC along the rungs. We 
looked for a magnetization plateau as follows. Motivated by the conditions in (|l|), we looked for a plateau 
at nis = 1/2 which would correspond to n = 1 in that equation, since S = 3/2. We also looked for 
plateaus at ms = and = 1, each of which would correspond to n = 2, i.e., a doubly degenerate state 
which has a period of two rungs. For a system with N sites, a given value of magnetization per rung, m^, 
corresponds to a sector with total equal to M = nisN/S. Using the infinite system algorithm, we found 
the lowest energies Eo{S^,N) in the three sectors = M + 1,M and M — 1. Then we examined the 
three plots of Eq/NJ versus 1/N and extrapolated the results up to the thermodynamic limit N oo. 
We fitted these plots with the formula Eq/NJ — Ci + Qi/N + bi/N'^, where the label i — 1, 2, 3 denotes 
the sectors M + 1, M and M — 1 . (We found that a quadratic fit in 1 /N matches the data much better 
than just a linear fit). In the thermodynamic limit, the values of the three intercepts should match 
since those are just the energy per site for the three states whose S'^'s differ by only 1. However, the three 
slopes ai are not equal in general. We now show that there is a magnetization plateau if ai -I- 03 — 2a2 has 
a nonzero value. Since the three energies Eq are computed without including the magnetic field term, 
the upper critical field hc+ where the states with Sz ~ M + \ and M become degenerate is given by 

h,+ {N) = Eo{M + l,N) - Eo{M,N). (4) 

Similarly, the lower critical field he- where the states with Sz — M and M — 1 become degenerate is 
given by 

K-{N) = Eo{M,N) - Eo{M-l,N). (5) 

We therefore have a finite interval Ah{N) — hc+{N) — hc-{N) in which the lowest energy state with 
5^ = M is the ground state of the system with N sites in the presence of a field h. If this interval has 
a nonzero limit as ^ 00, we have a magnetization plateau. Thus, in the thermodynamic limit, the 
plateau width Ah/ J is equal to oi -|- 03 — 2a2- 

We will now quote our numerical results for J/ J' — 1/3. For a rung magnetization of = 1/2, i.e., 
M — N/6, we found the three slopes to be equal to 3.77, —0.02 and —1.93; see Fig. 1. This gives the 
upper and lower critical fields to be 



= ai - 02 = 3.79 , = a2 - as = 1-91 

Ah hr+ — hr- 

— = + J — — = 1.88 . 



(6) 
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This is a sizeable plateau width, and it agrees with the exact diagonalization resultsEj and with the second- 
order perturbation expansion which will be discussed in the next section. For a rung magnetization of 
rUs = 1, we found the to be equal to 4.97, —0.24 and —5.43. Thus the upper and lower critical fields 
are 

Y = 5.21 , ^ = 5.19 , ^ = 0.02 . (7) 

Finally, for a rung magnetization of = 0, we need the energies of states with M = and M = ±1. 
Since the last two states must have the same energy, we have oi = 03 and it is sufficient to plot only 
£■0(0, A^) and Eo{l,N) versus 1/iV. We found ai and 02 to be equal to 0.39 and 0.34. This gives the 
upper and lower fields to be 

^ - 0.05, ^ = -0.05, ^ = 0.10. (8) 

The plateau widths given in and are rather small. In Fig. ^, we indicate the plateau widths 
Ah{N)/J as a function of 1/A^ for — 1/2, and 1. We will see that the LEH in the next section 
actually predicts that there should be no plateaus at = and 1. 

Next, we computed the various spin correlations for the 120-site system. We studied the spin densities 
{S^ „) where the chain index a = 1, 2, 3 and n is the rung index. [Due to the rotation invariance about 
the z axis, the other two spin densities (S'^„) must vanish]. For the plateau at rUs = 1/2, we found that 

(Sin) = (Sin) = 0.27 , (5i„) - - 0.04 , (9) 

for values of n in the middle of the system. The spin densities are shown in Fig. ^. 

We also examined several two-spin correlations which can be denoted by {S^ [S^ ^) and (<S'^;<S'f/„). 
For the zz correlations, it is convenient to subtract the product of the two separate spin densities; the 
subtracted zz correlations then go to zero for large rung separations |Z — n|, just like the H — correlations. 
At iris — 1/2, we found that all these correlations decay very rapidly to zero as \l — n\ grows. In fact, 
the fall offs were so fast that we were unable to compute sensible correlation lengths. All the correlation 
lengths are of the order of one or two rungs as can be seen in Figs. ^ and 

On the other hand, for the state at — 0, we found that all the two-spin correlations decay quite 
slowly. The decays are consistent with power law fall offs of the form A(— l)l'^"l/|/ — n\^. It is difficult to 
find 77 very accurately since the maximum value of |/ — n| is only 20; this is because we fixed one site to be 
in the middle of the chain (to minimize edge effects) , and the maximum chain length is 40 for our DMRG 
calculations. For rUs = 0, the exponent 77 for all the correlations was found to be around 1. There was 
no difference in the behaviors of the zz and H — correlations since this was an isotropic system; rUs — 
is the ground state if the magnetic field is zero. 

For the state at rUs = 1 (which is the ground state only for a substantial value of the magnetic field), 
we found that the H — correlations again decay quite slowly consistent with a power law. The exponents 
77 for the different H — correlations varied from 0.61 to 0.70 with an average value of 0.66; see Fig. ^ 
for an example. However, the zz correlations actually increased, rather than decreased, with increasing 
separation |/ — n\; see Fig. ^ We found that this is because of large edge effects. Since the magnetic 
field is particularly strong for the state with iris = 1, and sites at the ends have fewer neighbors coupled 
antiferromagnetically to them, they respond more strongly to the magnetic field than sites near the center 
of the system. This can be seen from Fig. ^ where the spin density „ shows a sharp increase towards 
the end of the chain (the rung index n is equal to 1 at the end) . 

We now summarize the properties of the three states studied with OBC along the rungs. The state 
with ms = 1/2 is characterized by a large gap to excited states and extremely short correlation lengths 
for spin correlations. The states at nis = and = 1 appear to have no gaps to excited states (within 
our numerical accuracy), and have slow fall offs of correlation functions consistent with power laws. 
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We now describe some low-temperature thermodynamic properties of the three-chain system obtained 
using DMRG. Ahhough DMRG is normally expected to be most accurate for targeting the lowest states 
in different sectors, earlier studies of mixjed spin chains have shown that DMRG is quite reliable for 
computing low-temperature properties alsocS. There are two reasons for this; the low- lying excited states 
generally have a large projection onto the space of DMRG states which contains the ground state, and 
the low-lying excitations in one sector are usually the lowest states in nearby sectors. 

We first checked that for systems with 12 sites, the results obtained using DMRG agree well with those 
obtained by exact diagonalization. We then used DMRG to study the magnetization, susceptibility and 
specific heat of 36-site systems with both OBC and PBC along the rungs. We first compute the partition 
function Z ~ exp [— /3(£'i — /i(S'^)i)] , where the sum is over all the states i in all the S'^ sectors, and 
/3 = l/ksT where kg is the Boltzmann constant. Then the magnetization is given by 

(M) = 1 y {Sn, e- /3 [ - is'h ] (10) 

i 

The magnetic susceptibility is related to the fluctuation in magnetization, 

X ^ Pi (M2) - (Af)2 ] , (11) 

and the specific heat is related to the fiuctuation in energy, 

^ = P'[ {E') - {Ef ] . (12) 

The plots of magnetization versus magnetic field for various temperatures are shown in Figs. ^ and |l^ 
for OBC and PBC, respectively, along the rungs. The temperature T is measured in units of J/kg- We 
see that the plateau at = 1/2 disappears quite rapidly as we increase the temperature. With OBC 
along rungs, the plateau has almost disappeared at T = 0.4 which is substantially lower than the width 
Ah/ J = 1.88. The plots of susceptibility in Fig. |ll|for OBC also show no surprises. The susceptibility is 
(exponentially) small at low temperatures in the region of the plateau because the magnetic excitations 
there are separated from the ground state by a gap. 

However, the specific heats shown in Figs. |l2| and |l^ demonstrate an interesting difference between 
OBC and PBC along the rungs. While it is very small at low temperatures for OBC, it is not small for 
PBC; further, it shows a plateau in the same range of magnetic fields as the magnetization itself. These 
two observations strongly suggest that the system with PBC along the rungs has nonmagnetic excitations 
which do not contribute to the magnetization or susceptibility, but do contribute to the specific heat. 
Figure ^ gives a more direct comparison between OBC and PBC along the rungs. The LEH of Sec. 
Ill will |djeaj"|ljz; show the origin of these excitations. Although these excitations were studied by previous 
author aljlljo, we believe that our specific heat plots prove their existence most physically. To show 
these excitations even more explicitly, we present in Fig. |l5| all the energy levels for a 12-site chain in 
the sector S*^ = 2 (i.e., = 1/2) using exact diagonalization. It is clear that the ground state is well 
separated from the excited states for OBC, but it is at the bottom of a band of excitations for PBC; 
these excitations are nonmagnetic since they have the same value of 5'^ as the ground state. 

We should point out that the rapid but small fluctuations seen in Figs. |ll| - in the susceptibility 
and specific heat at the lowest temperature of T = 0.1, are due to finite-size effects. Apart from a large 
plateau at uig = 1/2, a system with only 36 sites also has small plateaus for several values of m at zero 
temperature. These lead to small wiggles in the magnetization (M) at very low temperature. The wiggles 
get amplified in the susceptibility since it is equal to the first derivative, i.e., x — d{M)/dh. The specific 
heat shows low-temperature fluctuations for the same reason. 

We should mention here that a small plateau has been found at = for PBC along the rung£3'0. 
The half-width is given by hc+/J — 0.21 in the limit J'/ J oo. However, this plateau is not clearly 
visible in our low-temperature plots of magnetization and susceptibility. 
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III. LOW-ENERGY EFFECTIVE HAMILTONIANS 



A. General comments 



We will now discuss the LEH approach for studying the properties of spin ladders. There are two 
possible limits which may be considered. One could examine J' / J which corresponds to weakly 
interacting chains, and then direct| ly || U a e .techniques from bosonization and conformal field theory; this 
has been done in detail by othergHJ^Ollalla. We will therefore consider the strong-coupling limit J/ J' 
which corresponds to almost decoupled rungs. In that the LEH has been derived|-tO|^rst order in 

J/ J' for a three-chain ladder with PBC along the runggHS^Ell, and for a two-chain laddenla'Ea. 

We will derive the LEH for the three-chain model with OBC along the rungs and a two-chain model 
to second order in J/ J', and for the three-chain model with PBC along the rungs to first order. For the 
three-chain system with OBC and for the two-chain system, we find that the first-order LEH takes the 
form of the XX Z spin-1/2 madd-.in a magnetic field. A lot of information is available for this model 
through conformal field theoryll2l'E2l. In particular, the exponent rj for the correlation power laws can be 
read off from the first-order Hamiltonian. We will use the terms of second order in J/ J' only to determine 
the boundaries hc± of the various plateaus. The second-order terms should also give corrections to the 
exponent rj but we will not consider that problem here. For the three-chain model with PBC along the 
rungs, even the first-order LEH is sufhcicntly complicated that its properties are not well understood; 
however we will present the form of the LEH for completeness. 

We derive the LEH as follows. We first set the intrachain coupling J — and consider which of the 
states of a single rung are degenerate in energy in the presence of a magnetic field. In general, there will 
be several values of the field, denoted by ho, for which two or more of the rung states will be degenerate 
ground states. We will consider each such value of /iq in turn. The degenerate rung states will constitute 
our low-energy states. If the amount of degeneracy in each rung is d, the total number of low-energy 
states in a system with L rungs is given by L'^. (The number d depends both on the system and on the 
field hf). It is two for three chains with OBC along the rungs and for two chains, while it is three or four 
for three chains with PBC along the rungs. The form of the LEH depends crucially on this degeneracy). 
Next, we decompose the Hamiltonian of the total system as H ^ Hq + V, where Hq contains only the 
rung interaction J' and the field ho, and V contains the small interactions J and the residual magnetic 
field h — ho which are both assumed to be much smaller than J'. Let us now denote the degenerate and 
low-energy states of the system as pi and the high-energy states as qa- The low-energy states all have 
energy Eq, while the high-energy states have energies according to the exactly solvable Hamiltonian 
Hq- Then the first-order LEH is given, up to an overall constant, by degenerate perturbation theory, 



The calculation of the various matrix elements in Eqs. (|Tj) and can be simplified by using the 
symmetries of the perturbation V, e.g., translations and rotations about the z axis. 

Finally, if there is a state pi such that {pj\V\pi) = for all low-energy states j ^ i, then the unnormal- 
ized state pi is given, to first order, by 




(13) 



The second-order LEH is given by 




(14) 




(15) 



7 



This result will be used to compute the first-order changes in some quantities like the spin densities and 
correlations. 

Before ending this section, we would like to make a few_cnmments on the XX Z spin-1/2 model in a 
magnetic field since this will play an important role belowM^ES. Consider a spin-1/2 chain governed by 
the Hamiltonian 

H = ^ [ S^S^^i + S'^Sf^^i + A S^S^^i ] — h ^ . (16) 

n n 



where the anisotropy parameter A > — 1 . It is known that this system is gapped for h > 1 + A with all 
sites having — 1/2 in the ground state, and for ft, < —1 — A with all sites having Sz — ^1/2- For 
A < 1, these are the only two magnetization plateaus with m = ±1/2 per site. For A < 1 and h = 0, 
the two-spin correlations decay asymptotically as 



{Sq Sj^ ) 



4)" 



(-1)" 



\n\^/ri ' 



V = I: + - sin-i (A) . (17) 

Z TT 



On the other hand, for A > 1, there is an additional plateau at Wg = 0; there are two degenerate 
ground states which have a period of two sites consistent with the condition (|l|). Thus the invariance 
of the Hamiltonian under a translation by one site is spontaneously broken in the ground states. This 

is particularly obvious for A — > oo where the two ground states are H \ — • • • and — | h • ■ ■• The 

two-spin correlations decay exponentially for A > 1 and h = 0. 



B. Three-chain ladder with open boundary condition along the rungs 



We will decompose the Hamiltonian in (|^) as = Hq + V, where 

3 



a— 1,2 n a—1 n 

3 3 
= ^ E E Sa,n • Sa,„+1 " (^^ ^ /^o) E E ■ (18) 



a—1 n a—1 n 



We determine the field ho by considering the rung Hamiltonian Ho and identifying the values of the 
magnetic field ho where two or more of the rung states become degenerate. 

The eight states in each rung are described by specifying the 5^ components (+ and — denoting -t-1/2 
and —1/2 respectively) of the sites belonging to chains 1, 2 and 3. For instance, the four states with 
total S — 3/2 are denoted by |1), |4), where |1) = | + ++) and the other three states can be obtained 
by acting on it successively with the operator — J^a'^a- These four states have the energy J'/2 
in the absence of a magnetic field. There is one doublet of states |5) and |6) with 5* = 1/2, where 

|5) = [ 2 I H h) - I - ++) - I + +-) ]/\/6 and |6) S^\5). These have energy - J'. Finally, there is 

another doublet of states |7) = [ | -I- H — ) - | - ++) ]/V2 and |8) S~\7) which have zero energy. It is 
now evident that the state |1) with — 3/2 and the state |5) with — 1/2 become degenerate at a 
magnetic field Hq — 3J'/2, while states |5) and |6) are trivially degenerate for the field Hq — 0. We will 
now examine these two cases separately. 
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For ho = 3J'/2, the low-energy states in each rung are given by |1) and |5), while the other six are 
high-energy states. We thus have an effective spin-1/2 object on each rung n. We may introduce three 
spin-1/2 operators (5*^, S"^, S^) for each rung such that = S^± iS^ and have the following actions: 

S+ |1)„ = , S+ |5)„ = |1)„ , 

■^n 11)" ~ ' 1^)" ~ ^ J 

|l>n - \ |1>» , |5)„ = - i |5)„ . (19) 



Note that the state which has a |1) on every rung, i.e., |111 • • •), is just the state with rung magnetization 
rris = 3/2 corresponding to the saturation plateau. The state with a |5) on every rung corresponds to the 
rris = 1/2 magnetization plateau. The LEH we are trying to derive will therefore describe the transition 
between these two plateaus. 

We now turn on the perturbation V in (|8|) with the assumption that J and h — ho are both much 
smaller than J'. We can write V = J2n ^n,n+i, where 

3 1 ^ 

Vn,n+1 = J E Sa.„ •Sa,„+i - -{h - ho) [ S^^ + S^^^, ] . (20) 

a— 1 a— 1 



The action of Vn,n+i on the four low-energy states invol ving rungs n and n -\- 1 can be obtained after a 
long but straightforward calculation. We then use Eq. ( [l4| ) and find that the LEH to second order in 
J/ J' is given, up to a constant, by 

-^e// = J E [ 5'^5'^_|_i -I- S'^S'^^j^i ^2 ^ 72,]'^ ^n^n+l ] 

n 

5J2 



^ 1 



18J' 

n 

3J' J 29. 
1 2 ^ 72J 



(/^-^-^- t^)E^n. (21) 



where we have substituted ho — 3J'/2. Note that the terms of order J only involve two neighboring sites. 
The LEH up to that order simply describes an XX Z model with anisotropy A = 1/2 in a magnetic field 
h — 3J'/2 — J/2 [see (p^j. Some of the terms of order J^/ J' involve three neighboring sites; this makes 
the model unsolvable by the Bethe ansatz at this order. 

We will now use ( ^l|) to compute the values of the fields hi and /12 where the states with all rungs 
equal to |1) and all rungs equal to |5) respectively become the ground states. We can then identify hi 
with the lower critical field he- for the plateau at = 3/2, and with the upper critical field hc+ for 
the plateau at to^ = 1/2. [Recall the definition of upper and lower critical fields around Eqs. (0) and 




To compute the field hi, we compare the energy Eo of the state with all rungs equal to |1) with the 
minimum energy Emin{k) of a spin- wave state in which one rung is equal to |5) and all the other rungs 
are equal to |1). A spin wave with momentum k is given by 

= 4= E e''" : (22) 
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where |5„) denotes a state where only the rung n is equal to |5). The spin-wave dispersion, i.e., Lu{k) = 
E{k) — Eq, is found from to be 

.,, ^ . , 1 29 J , , , 3 J J 29 , . , 

u;{k) = J ( cosfc - - + —^) + {h- — - -- -j^). (23) 

This is minimum at fc = tt and it turns negative there for h < hi, where 

hi ^ ^- + 2J. (24) 

This is therefore the transition point between the ferromagnetic state |111 • • •) and a spin- wave band lying 
immediately below it in energy. 

Similarly, we compute the field /12 by comparing the energy Eq of the state with all rungs equal to |5) 
with the minimum energy Eminik) of a spin wave in which a |5) at one rung is replaced by a |1). For a 
spin wave with momentum k, the dispersion cj(fc) = E{k) — Eq is found to be 

1 29J , J2 2 5 , , , 3J' J 29J2 , ,^ , 

= J (cask - 2 + ^ ) + 77 ( 9 - cos2fc) 2 - ^ ) . (25) 

This is minimum at fc = tt and it turns positive there for /i > /i2, where 

3 /' 3 P 

This marks the transition between the state |555- • •) and the-spin-wave band. Equation ( p6| ) agrees to 
this order with the higher-order series given in the literature^. Note that the second-order result gives 
h2/J = 3.75 for J/ J' = 1/3, compared to our DMRG value of K+/J = 3.79 in (|). 

From the first-order terms in (^), we can deduce the asymptotic form of the two-spin correlations. 
From (|l^, we see that the exponent 77 — 2/3 for A = 1/2. Although this is the exponent for the H — 
correlation of the effective spin-1/2 defined on each rung, we would expect the same exponent to appear in 
all the correlations {S^iSf^^) studied by DMRG in the previous section, regardless of how we choose the 
chain indices a,b — 1, 2, 3. We now see that the analytically predicted exponent of 2/3 agrees quite well 
with the numerically obtained exponents which lie in the range 0.61 to 0.70. Incidentally, this agreement 
implies that J/ J' = 1/3 is sufficiently small so that the second-order terms do not significantly affect the 
correlation exponent. 

Finally, we can use the first-order wave function given in (^) to compute the spin densities and short- 
distance two-spin correlations. As examples, we quote the results for spin densities and some of the 
nearest-neighbor spin correlations for the plateau at ms = 1/2. We will give the first-order expressions 
and their values for J/ J' = 1/3, followed by the numerical values obtained by DMRG. 

1 4J 

{SU = {SU = 3 - ^ = 0.28 vs. 0.27 from DMRG , 

1 8(7 

(Sin) - - 6 + 277 = " ~ ' 

(SlnSln+i) = \ - = 0-07 vs. 0.07 , 

1 4 / 

{SlnSln+i) ^ ^ - ^ = - 0-02 vs. - 0.02 , 

{SLSi,n+i) - - {j, - - 0-06 vs. - 0.08 , 

2 7 
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We will now consider the LEH at the other magnetic field ho = where the rung states |5) and |6) are 
degenerate. We take these as the low-energy states and introduce new efi^ective spin-1/2 operators for 
each rung with actions similar to Eqs. ([T9|), except that we replace |1) and |5) in those equations by |5) 
and |6). We again compute the action of the perturbation V on the low-energy states. We then deduce 
the second-order LEH to be 

^e// ^ X! [ ( ^ 9j7 ) ^" ■ ^"+1 ^ 27J' ^" ' ^"^^ ] ^ ^ ■ (28) 

n n 

This Hamiltonian describes the transition between the magnetization plateaus at = 1/2 and = 
— 1/2; since these plateaus are reflections of each other about zero magnetic field, it is sufficient to study 
one of them. By a calculation similar to the one used to derive (^4|), the field hi can be found from the 
dispersion of a spin wave in which one rung is equal to |6) and all the other rungs are equal to |5). The 
dispersion is 

j2 g j2 

u;{k) = h + [J - — ) { cosfc - 1 ) + ^ ( 1 _ cos2fc ) . (29) 



This gives 



2/2 

hi = 2J - —. (30) 



This is the lower critical field hc~ of the = 1/2 plateau. For J/ J' — 1/3, the second-order result 
gives 1.93 versus the DMRG value of 1.91 in (^. The Hamiltonian ( p8| ) describes an isotropic spin-1/2 
antiferromagnet with a weak ferromagnetic next-nearest-neighbor interaction. From the comments at the 
end of the previous section, we see that this model only has the two saturation plateaus at = ±1/2, 
and no other plateau in between. For /i = 0, the two-spin correlations decay as power laws with the 
exponent rj = 1 [see (p7|)]. 



C. Three-chain ladder with periodic boundary condition along the rungs 

In this section, we will present the first-order LEH for the Hamiltonian (H) with PBC along the rungs. 
The LEH will turn out to be somewhat complicated. We will not study their properties in any detail, 
but will limit ourselves to a few comments. As in the case with OBC along the rungs, there are two 
different LEH to be considered here because there are two values of the magnetic field where there are 
degeneracies. We again begin with a description of the eight states on each rung. The four states with 
S = 3/2 are the same as the states |1), |4) introduced in the previous section, except that they now 
have energy 3J'/4 in the absence of a field. The doublet states have to be chosen differently now in 
order that they be eigenstates of the periodic rung Hamiltonian. We choose two of the doublet states 
to be |5') = [\+ +-) \+ -+) +UJ \- ++) ]/^/3 and |6') - S'-|5'), where uj = exp(i27r/3). 

These two states have momenta 27r/3 along the rung (right moving). The other two doublet states are 

|7') = [ I + +-) h) + I - ++) ]/V3 and |8') - S-\l') with momenta -27r/3 (left moving). 

All these four states have energy — 3J'/4. This extra degeneracy (which is twice the degeneracy of the 
doublets for OBC along the rungs) leads to a more complicated LEH as we will see. 

We now note that for a field /iq = 3J'/2, the three states |1), |5') and \T) become degenerate. We 
now introduce seven operators R^,L^,t^ and for each rung with the following nonzero actions on 
the three low-energy states, 

Rn |5')« = |l)n , R- \l)n = |5')„ , 
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L+ |7')„ = |1)„ , L- |1)„ = |7')„ , 
T+ |7')„ = |5')„ , T- |5')„ = |7')„ , 

< |1)„ = |1)„ , < |5')„ - - |5')„ , < |7')„ = - \l')n ■ (31) 



All the actions (of operators on states) not mentioned in Eqs. ( |3l| ) are assumed to give zero. We thus 
observe that there are five magnetic operators L^, and which change or measure the of a state, 
and two nonmagnetic operators which do not change 5*^ but simply interchange the right and left 
moving states. 

We then find that the first-order LEH is given, up to a constant, by 

n 

+ ^2 [ '^n'''n+l + '''n '''n+1 ] + ^ ^-n^n+l 

n n 

1 / T 3t/' 2t/ , X — ^ ^ /^^\ 

n 

We can now find the magnetic field hi at which the the ferromagnetic state |111 • • •) crosses over to the 
minimum energy spin- wave state (in which a |1) is replaced by either a |5') or a |7') on exactly one rung). 
The spin-wave dispersion is 

3 J' 

(jj(fc) — h — — J + J cosfc . (33) 



We thus see that hi — 3J'/2 + 2 J just as for OBC along the rungs. Thus the lower critical field he- of 
the saturation plateau — 3/2 has the same value for OBC and PBC along the rungs. 

Below some field /12 (which seems rather hard to find analytically), the low-energy eigenstates of (^) 
will not have the state |1) on any rung; only the states |5') and |7') will appear. This gives us the 
magnetization plateau vris = 1/2. However, this plateau has a large number of nonmagnetic excitations 
described by the Hamiltonian 

Heff = ^ Yl [ + ] (34) 

n 

which may be obtained from ( |32| ) by omitting the state |1) on all the rungs. Equation (^J) has the form 
of ( p^ ) with A = 0, and is therefore exactly solvable; at low temperature, it has a specific heat which 
grows linearly with T . The situation is therefore quite different from the case of OBC along the rungs 
where the nig = 1/2 plateau consists of a single state in which every rung is in the state |5); all other 
states are separated by a gap, hence the specific heat goes to zero exponentially at low temperature. 

Finally, we examine the LEH at the field h^ — Q-wheie the four doublet states |5'), • • • , |8') become 
degenerate. This has been discussed in detail earlierU'ElEj. On each rung, we introduce effective spin-1/2 
operators which change or measure 5^, and the two nonmagnetic operators which interchange the left 
and right moving states. Then the LEH is 

^e// = I E [ 1 + 4 ( r+r-^i + r-T+^i ) ] S„ • S„+i ~ hY^Sl. (35) 

n n 

This also appears to be nonexactly solvable but it has been studied numerically00. It has a small 
plateau at rus = where there are two degenerate ground states, each with a period of two rungs. Above 



12 



some magnetic field hi (whicli is again hard to calculate analytically from (pSj), this model crosses over 
to the rUg — 1/2 plateau where the rungs can only be in the two = 1/2 states |5') and |7'). We see 
from (^5|) that these two states are again governed by the Hamiltonian in (^). 

The phase diagrams of the three-chain ladder for OBC and PBC along the rungs are shown as functions 
of J'/J and h/J in Figs. 6 (b) and (c) in Ref. [11]. We observe that the plateaus with nis = 1/2 and 
■nis —3/2 (called M = 1/3 and M = 1, respectively, in Ref. [11]) have large regions of stability for both 
OBC and PBC. The plateau with iris — exists only in the case of PBC along the rungs, and it has a 
small region of stability close to h/ J — 0. 



D. A two-chain ladder 

In this section, we will use the LEH approach to study a two-chain spin-1/2 ladder with the following 
Hamiltonian, 

2 

H = J' Si^n • S2,n + J2 Sa,n ' Sa^„+i 

n a— 1 n 

2 

-I- 2 Ji Si^„ • S2,„+i - E E ^-'.n ' (36) 

n a—1 n 



as shown in Fig. |T^. The model may be viewed as a single chain with an alternation in nearest-neighbor 
couplings J' and 2Ji (dimerizatinn)., and a next-nearest-neighbor coupling J2 (frustration). Equation 
( p^ ) has been studied extensivelyEIIo. More recently, it has been studied from the paiHt-pf view of 
magnetization plateaus using a first-order LEH, bosonization and exact diagonalizationEjEaEa . We will 
therefore limit ourselves to deriving the second-order LEH and making a few other comments. 

We begin by setting Ji = J2 — 0, and studying the four states on each rung. These are specified by 
giving the configurations ± of the spins on chains 1 and 2 as follows. The three triplet states with S' = 1 
are denoted as |1), |2) and |3), where |1) = | + +) and the other two states are obtained by acting on it 
successively with . These three states have energy J'/4 in the absence of a magnetic field. The singlet 
state \A) = [\ + -)-\- +)]/V2 has energy — 3J'/4. The states |1) and |4) become degenerate at a field 
ho = J' . We now develop perturbation theory by assuming that Ji, J2 and h — ho are all much less than 
J' . The perturbation is V = J2n ^n,n+i where 

2 

Vn.n+l = J2 ^a,n ' Sa,n+1 + 2Ji Si^„ • S2,„+l 

0=1 

1 ^ 

-^ih-ho)J2[Sl„ + Sl„+i]. (37) 



The actions of this operator on the four low-energy states of a pair of neighboring rungs can be easily 
obtained. We now introduce effective spin-1/2 operators S„ on each rung which act on the two low-energy 
states. The second-order LEH is then found to be 

J2 

n 

, 1 / T , T , 2 3(Ji - J2)^ Sr qzqz 
+ 2^ + Jl + —, j 2^ 

n 
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72-1 1 
2JI 2-^ '-'njWn-l'-'n+l i '^n-l'-'n+lJ + I r, + "^n-l JWn'-'n+l + '-'n'-'n+l) 



4J' 



We now compute the field hi above which the state |111 • • •) becomes the ground state. The dispersion 
of a spin wave, in which one rung is equal to |4) and all the others are equal to |1), is given by 

J2 J2 

cj(fc) = - J' - Ji - J2 - + (-^2 - Ji) cos/c + ^ cos2fc . (39) 

By minimizing this as a function of k in various regions in the parameter space (Ji, J2), and then setting 
that minimum value equal to zero, we find that hi is given by 

2 P 

hi ^ J' + 2Ji if J2 < Ji - , 

-J + Ji + J2 + jj- + ^ It ^1 - — < ^2 < Ji + — , 

2 P 

= J' + 2J2 if J2 > Ji + -yr • (40) 

This is the lower critical field he- of the saturation plateau with magnetization rUs = 1 per rung. Similarly, 
we can find the field /12 from the dispersion of a spin wave in which one rung is equal to |1) and the rest 
are equal to |4). The dispersion is given by 

u{k) = - h + J' + - ^ + (J2-Ji-^) cosfc - ^^'-f^" cos2fc. (41) 

By setting the minimum of this equal to zero, we find that /12 is given by 

h, = f + ^-^^^j^ + J2 - Ji - ^ if J2 < Ji + ^ , 

= J' + 2/'^' - J2 + Ji if J2 > Ji + ^ . (42) 

This is the upper critical field hc+ of the saturation plateau with magnetization rUs = per rung. 

Finally, we can see that the first-order terms in (38) are of the same form as the XX Z model in (p^). 



We can always make the coefficient of the first term in (^8|) positive, if necessary by performing a rotation 
Sn ^ i-'^Y'S^, SI {-ITSI and S^. We then get a first-order Hamihonian of the form 

n n 

-{h- J' -'^ -'^)Y^S:. (43) 
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This is an XX Z model with 



From the comments in Sec. Ill A, we therefore see that the two-chain ladder will have an additional 
plateau at = 1/2 for A > 1, i,«,, if J2 + Ji > 2| J2 — Ji|. In particular, A = 00 for J2 = Ji (this is 
called the Shastry-Sutherland lineEJ); the = 1/2 plateau will then stretch all the way from the upper 
critical field of the — plateau to the lower critical field of the — 1 plateau. This can be seen 
in Fig. [ij which is taken from Ref. [16]; the dimerization parameter a in that figure is related to our 
couplings by J' = 1 + a and 2Ji = 1 — a. Note that the nis = 1/2 plateau is particularly broad at 
a = 0.6, i.e., J2 = Ji = 0.2, and that it actually touches the nis = 1 plateau on the right. The fact that 
it does not extend all the way up to the rus = plateau on the left is probably because we have ignored 
the second-order terms in (p8|) which lead to deviations from the XX Z model. 



IV. SUMMARY AND OUTLOOK 



We studied a three-chain spin-1/2 ladder with a large ratio of interchain coupling to intrachain coupling 
using the DMRG method and a LEH approach. For both OBC and PBC along the rungs, we found a 
wide plateau with rung magnetization given by = 1/2. For the case of OBC, the two-spin correlations 
are extremely short-ranged, and the magnetic susceptibility and specific heat are very small at low 
temperature in the plateau. All these are consistent with the large magnetic gap. At other values of 
m, the two-spin correlations fall off as power laws; the exponents can be found by using the first-order 
LEH which takes the form of an XX Z model in a longitudinal magnetic field. For the case of PBC, the 
magnetic susceptibility is again very small at low temperature in the plateau. However the specific heat 
goes to zero much more slowly which dramatically shows the presence of nonmagnetic excitations. This 
can be understood from the LEH in (|3^) which is an XY model. Finally, we used the LEH approach to 
study a two-chain ladder with an additional diagonal interaction. In addition to a plateau at nis = 0, this 
system also has a plateau at ms = 1/2 for certain regions in parameter space. The rris = 1/2 plateau is 
interesting because it corresponds to degenerate ground states which spontaneously break the translation 
invariance of the Hamiltonian. This can be understood from the LEH which, at first-order, is an XX Z 
model with A > 1. 

An interesting problem for the future may be to take the second-order terms in the LEH presented 
in Sees. HI B and HI D, and to compute the corrections produced by them in the exponents of the 
correlation power laws. This would require us to study the effects of a perturbation to the XX Z spin-1/2 
chain. This may not be difficult to do analytically since the XX Z model is integrable and exactly solvable 
by the Bethe ansatz. 

The quantization condition for magnetization given in (|l]) is reminiscent o£,the quantum Hall effect 
where the Hall conductivity shows plateaus as a function of the magnetic fielcE2l. However, it is not clear 
if the magnetization quantization is as insensitive to disorder as the conductivity quantization is known 
to be. Although a magnetization plateau may be expected to survive small amounts of disorder (e.g., if 
the disorder strength is much smaller than the energy gap), there seems to be no fundamental physical 
principle, analogous to gauge invariance in the quantum Hall system, why the value of the magnetization 
should remain fixed at a simple rational value. In fact, the derivation of {M assumes translation invariance 
of the Hamiltonian which is certainly broken by disorder. It would therefore be interesting to study this 
issue, for instance, by allowing a small amount of disorder in the couplings of the spin ladder models 
discussed in this paper. 
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Figure Captions 

1. The energy/site in units of J vs l/N at the = 1/2 plateau, for J/ J' = 1/3. The curves indicate 
quadratic fits for (a) Eo{M + 1,7V), (b) Eo{M,N), and (c) Eo{M - 1,7V). 

2. Plateau widths vs 1/7V for (a) = 1/2, (b) rus = 0, and (c) rUg = 1. 

3. Spin densities at the = 1/2 plateau for J /J' = 1/3. The upper points (circles) denote the top 
chain a = 1, while the lower points (triangles) denote the middle chain a = 2. n = 1 and 20 denote the 
end and middle rungs respectively. 

4. Correlation function < 5'^;<S'^„ > at the rUs = 1/2 plateau for J/ J' = 1/3. 

5. Correlation function < iSf „_ > at the = 1/2 plateau. 

6. Correlation functions < 5'^;S'^„ > in the = 1 state for J/ J' = 1/3. 

7. Correlation function < S'|;S'|„ > in the m,, = 1 state. 

8. Spin densities in the mg = 1 state for J/ J' — 1/3. The upper points (circles) denote the top chain 
a = 1, while the lower points (triangles) denote the middle chain o = 2. n = 1 and 20 denote the end 
and middle rungs respectively. 

9. Magnetization vs magnetic field for 36 sites, with OBC along rungs for J/ J' = 1/3. 

10. Magnetization vs magnetic field for 36 sites, with PBC along rungs for J/ J' = 1/3. 

11. Susceptibility vs magnetic field for 36 sites, with OBC along nmgs. 

12. Specific heat in units of ks vs magnetic field for 36 sites, with OBC along rungs for J/ J' = 1/3. 

13. Specific heat in units of ks for 36 sites, with PBC along rungs. 

14. Comparisons of specific heat and susceptibility of the 36-site systems with OBC and PBC along the 

rungs. 

15. Comparison of the energy spectra in units of J of the 12-site system with OBC and PBC along the 
rungs. The energies in the Sz = 2 sector are shown for J/ J' = 1/3. 

16. Schematic diagram of the two-chain ladder with an additional diagonal interaction. The labels 1 and 

2 denote sites in the upper and lower chains respectively. 

17. Phase diagram of the two-chain ladder as a function of h and a for J2 = 0.2. In our notation, 
J' = 1 + a and 2Ji = 1 — a. The numbers 0, 1/2 and 1 in the figure correspond to the values of rUs at 
the plateaus. Reproduced with permission from Ref. [16]. 
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